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ABSTRACT 

Using the Oxford Short Wavelength Integral Field specTrograph, we investigate radial 
variations of several initial mass function (IMF) dependent absorption features in M31 
and M32. We obtain high signal-to-noise spectra at six pointings along the major axis 
of M31 out to ~ 700" (2.7 kpc) and a single pointing of the central 10 pc for M32. 
In M31 the sodium Nal A8190 index shows a flat equivalent width profile at ~ 0.4 A 
through the majority of the bulge, with a strong gradient up to 0.8 A in the central 
10"(38pc); the Wing-Ford FeH A9916 index is measured to be constant at 0.4 A for 
all radii; and calcium triplet CaT A8498, 8542, 8662 shows a gradual increase through 
the bulge towards the centre. M32 displays flat profiles for all three indices, with FeH 
at ~ 0.5 A, very high CaT at ~ 0.8 A and low Nal at ~ 0.1 A. We analyse these data 
using stellar population models. We find that M31 is well described on all scales by 
a Chabrier IMF, with a gradient in sodium enhancement of [Na/Fe] ~ +0.3 dex in 
the outer bulge, rising within the central 10" to perhaps [Na/Fe] ^ +1.0 dex in the 
nuclear region. We find M32 is described by a Chabrier IMF and young stellar age in 
line with other studies. Models show that CaT is much more sensitive to metallicity 
and [a'/Fe] than to IMF. We note that the centres of M31 and M32 have very high 
stellar densities and yet we measure Chabrier IMFs in these regions. 
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1 INTRODUCTION 

One of the most rigorous methods of constraining galaxy 
formation and evolution is through analysis of their stellar 
spectra. Stars are the directly observable baryonic compo¬ 
nent of galaxies, which habour information about their for¬ 
mation. Integrated light from the unresolved stellar compo¬ 
nent of galaxies thus presents a direct probe of their compo¬ 
sition. The initial mass function (IMF) is one of the most im¬ 
portant parameters to constrain our understanding of galaxy 
formation. This function sets the mass distribution in stars 
and affects the total mass-to-light ratio of the galaxy. It has 
been one of the most difficult functions to constrain for ex¬ 
ternal galaxies, as direct star counts are not possible. The 
IMF is defined as the number density of stars n per unit 
mass m formed at one time. The work of Salpeter (1955) 
(and later Scalo 1986) showed that the IMF £(ra) takes a 
simple functional form in the Milky Way 


£(ra) — — km 

dm 
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and the exponent x takes a value of x — 2.35, down to masses 
of ^ 1 M 0 . The low mass end of the IMF is especially chal¬ 
lenging to constrain from integrated light, as stars of mass 
< 0.5 M 0 contribute only ~ 1 per cent of the bolometric 
luminosity of the galaxy. However, several studies extended 
the form of the IMF below 1 M 0 , and suggested either a flat¬ 
tening of the slope ( Miller &; Scalo|1979 ) or even a turnover 
around 1 M 0 ( Kroupa|2001 Chabrier|2003 ). 

It has long been assumed that the form of the IMF 
is universal for all galaxy types, masses and ages. Indeed, 
most galaxy formation and large scale structure simulations 
simply invoke a single form of the IMF for all simulations. 
Throughout the second half of the 20th century, several pa¬ 
pers highlighted and detected specific stellar absorption fea¬ 
tures that are prominent in the spectra of low mass stars, but 
absent in high mass counterparts (e.g. Wing-Ford band FeH 


A9916: Wing & Ford 1969[|Schiavon, Barbuy &; Singh|1997[ 
Cushing et al.||2003 sodium Nal A8190: |Faber &; French] 


1980 


Schiavon et al. 1997). A lack of both high signal- 


to-noise and detailed stellar populations models, coupled 
with the challenges of spectroscopy in the far red and near- 
infrared, meant that arguments for a varying IMF were in¬ 
conclusive. However, recent progress by Ivan Dokkum & Con- 


roy (2010) and Conroy & van Dokkum | ( |2012 b| has shown 
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that it is possible to detect and model these absorption fea¬ 
tures in integrated light from massive early type galaxies 
(ETGs) with very high signal-to-noise spectra. They de¬ 
tected the Nal doublet feature at 8190 A, the calcium triplet 
CaT around 8500 A, and the FeH Wing-Ford band at 9916 A 
in spectra from the centres of 38 giant ETGs. Their results 
have shown that these galaxies are well modelled by bottom- 
heavy (x > 2.35) IMFs and that the IMF becomes more 
bottom-heavy with increased velocity dispersion. The com¬ 
plimentary approaches of lensing/dynamical modelling to 
constrain the form of the IMF in ETGs (e.g. strong lensing; 
|Treu et al.|2 010 ; Smit h, Lucey & Conroy]2015] and integral 


field kinematics; Cappellari et al. 2012), has also provided 
consistent constraints with spectroscopic methods (although 
quantitive differences exist; see Smith 2014). Further, work 
by Geha et al. (2013) using HST deep imaging has shown ev¬ 
idence for a very bottom-light IMF in very diffuse galaxies, 
lending further weight to non-universality of the IMF. 

While the work by |Conroy fe van Dokkum| and others 
shows a paradigm shift is required from the assumption of a 
constant IMF for all galaxies, their results are only represen¬ 
tative of the central regions of each galaxy. Searching for ra¬ 
dial gradients in these absorption features and the IMF will 
provide more conclusive results for use in galaxy modelling 
and stellar population analysis. Recent results by |Martfn-| 


Navarro et al. (2015) have indicated a gradient to the IMF 


slope in a massive early type galaxy with large central ve¬ 
locity dispersion. However, work by Smith & Lucey (2013) 
indicated the presence of massive ETGs with bottom-light 
IMFs. These results show that our knowledge of the IMF is 
still uncertain. 

While all of the recent IMF work has looked at ETGs, 
it would be interesting to compare to the bulges of spi¬ 
ral galaxies. M31 is an ideal target for high S/N resolved 
spectroscopy due to its proximity. It is an archetypal spiral 
galaxy and the central bulge has had its stellar populations 


extensively studied, both using spectroscopy (e.g. Davidge 
Saglia et al.|2010 ) and photometry (e.g. Stephens et al. 


1997 


2003 


Sarajedini &; Jablonka 2005). Determining the pres¬ 


ence of IMF gradients will help provide a more comprehen¬ 
sive understanding of spiral bulges and the process behind 
IMF variations. 

The form of the IMF in the nuclear region of M31 
has been the subject of much debate throughout the sec¬ 


ond half of the 20th century (see Conroy & van Dokkum 
2012b] section 4.2). The core has been known to have promi¬ 


nent sodium absorption and this has been interpreted as an 


increased dwarf abundance by some (e.g. Faber & French 


1980) and sodium enhancement by others (e.g. Alloin & Bica 


1989). 


In this paper we present three studies on: 

a) the variations of far red absorption features in M31 
through the bulge out to the old disc at ~ 700" (2.7kpc); 

b) the variations of these features through the inner bulge 
to the nuclear region of M31 where there has been much 
historic debate; 

c) these features in the central 8" (30 pc) of M32. 

We use the definitions from Kormendy & Bender (1999) 
and Tempel et al.| ( |20lT] ) to describe the components of M31. 
The nuclear region is defined as r < 4"(15pc), the bulge as 
4 " < r < 450"(1.7kpc), and the disc as r > 450". We also 
refer to the region between 4" < r < 10" as the inner bulge. 


Table 1 . Observational fields of M31 and M32. Observations of 
M31 are labelled as M31_l through to M31_6 in increasing radii 
from the centre. We also give the positions of the sky fields. 


Field 

RA 

DEC 

Date 

M31_l 

00:42:44.57 

+41:16:05.7 

03-12-2011 

M31_2 

00:42:46.82 

+41:16:41.3 

23-09-2012 

M31_3 

00:42:51.04 

+41:17:38.3 

23-09-2012 

M31_4 

00:42:58.19 

+41:19:17.8 

02-10-2014 

M31_5 

00:43:10.80 

+41:22:03.7 

23-09-2012 

M31_6 

00:43:28.30 

+41:25:52.8 

23-09-2012 

M32 

00:42:41.86 

+40:51:57.0 

23-09-2012 

SKY_M31 

00:44:12.04 

+40:48:05.6 


SKY_M32 

00:46:15.27 

+41:01:28.3 



We analyse our results using the latest stellar popula¬ 
tion synthesis models from Vazdekis et al. (2012 hereafter 
V12) and Conroy &; van Dokkum (2012a hereafter CvD12). 
We further investigate a method of disentangling the degen¬ 
eracies in age and metallicity by using colour as a proxy for 
age, in the absence of the standard H/3 absorption line. 

This paper is organised as follows: section 2 details our 
observations and data reduction procedures. In section 3 
we present the absorption feature strengths, and section 4 
contains our stellar population analysis using simple stellar 
population (SSP) models and colour-index maps. We dis¬ 
cuss our results in section 5 and conclude our work in sec¬ 
tion 6. In this paper we adopt a distance to M31 of 785 kpc 
as measured by McConnachie et al. (2005), which gives a 
scale of 3.8pcarcsec _i . We further adopt the IMF naming 
convention used by CvD12, namely a Chabrier IMF corre¬ 
sponding to that defined in Chabrier (2003) for the disc of 
the Milky Way; a Salpeter IMF defined as x — 2.3 for V12 
and x = 2.35 for CvD12; and a bottom-light IMF as defined 
(2008). 


van Dokkum 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Overview 


Observations of M31 and M32 were obtained over the nights 
of 3rd and 5th December 2011; 23rd and 24th September 
2012; and 2nd October 2014 using the Short Wavelength 


Integral Field specTrograph (SWIFT) instrument (Thatte 


et al. 2006) on the Palomar 200 inch (5.1 m) telescope. 
Fig. [l shows the position of each pointing. Observations 
were taken at the 235 mas spaxel -1 spatial scale covering 
a field of 10" X 21". The spectra cover the wavelength range 
6300 — 10400 A with a dispersion of ~ 2 A FWHM and a 
sampling of 1 Apix -1 . Table [l] lists the pointings and dates 
of the observations. The seeing was around 1.5" for all ob¬ 
servations. 

The data were reduced using the SWIFT data reduction 
pipeline, written in iraf (Houghton, in preparation). The 
pipeline handles all the standard reduction processes of bias 
subtraction, flat fielding, wavelength calibration, error prop¬ 
agation as well as IFS specific features of illumination correc¬ 
tion and cube reconstruction. Cosmic rays were detected and 
removed using the LACOSMIC routine (van Dokkum|2001). 

We obtained Sloan Digital Sky Survey (SDSS) g — r 
colour data from Saglia et al. (2010) and | Temp el et al. 
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Figure 1 . STScI/Digital Sky Survey (DSS) image of M31 and M32 showing the locations of the SWIFT fields (blue boxes). The M31 
fields are named in order of increasing radius from M31_l through to M31_6. Note that the boxes only show the locations of the fields, 
and are not representative of the size of the SWIFT field of view. 


(2011). Saglia et al. 


for which we used 


covered all M31 pointings except M31_6 
Tempel et al. (2011). Sodium D (NaD 


A5893) index measurements for M31 and M32 were obtained 
from |Davidge (1991 1997). We applied corrections to the 
Davidge NaD data based upon the CvD12 models to ac¬ 
count for the measurements being at a lower resolution than 
a = 200 km s -1 (see section 2.5). 


2.2 Sky subtraction and spectra extraction 


Sky subtraction was performed to first order by subtracting 
sky frames observed adjacent in time to each science frame. 
We were careful to choose sky fields sufficiently far away as 
to avoid contamination from M31. The positions are given 
in table ^ We extracted spectra in two ways. 

For fields M31_2 to M31_6 we summed up all spaxels 
within the field to obtain a median spectrum of maximal 
S/N; in this case we first normalised by the median image 
to remove the increased weighting of brighter regions over 
fainter ones. 

For M31_1 and M32 fields we performed sector bin¬ 
ning ( |Houghton et al.||201 3) to obtain a 2D map of spectra 
with S/N > 100 pixel - . This gave ~ 200 individual spectra 
within a radius of 10 x/ for M31_l and ~ 170 spectra within 
8" for M32. The numbers are slightly lower for the FeH op¬ 
timised spectra due to more demanding S/N requirements 
for that region of the spectrum. 


2.3 Telluric correction 

Atmospheric telluric absorption is prominent in both the 
Nal and FeH regions. We observed A0V stars during the 
observations to act as telluric standards. These stars have 
strong Paschen lines around 8500 A which we remove using 
a dedicated routine written in idl. The telluric spectra are 
divided out by model A0V spectra provided by R. Kurucz, 
using a version of the amoeba algorithm allowing for velocity 
shifts and stellar rotation. Our corrected telluric spectra are 


accurate to around the 0.5—1 per cent level (see Fig. Al 
appendix). 


2.4 Kinematic fitting and second-order sky 
subtraction 


We fit the kinematics of the spectra with the penalised 


pixel fitting routine (ppxf, Cappellari & Emsellem 2004) 
using, as template spectra, the stellar population models 
from CvD12. Due to the presence of skyline residuals in the 
spectra, we found it necessary to perform a simultaneous 
second-order sky subtraction while extracting the kinemat¬ 
ics ( Weijmans et al.|2009 ). This was accomplished using the 
techniques described in Davies (2007): rather than scaling a 


single sky spectrum covering multiple OH vibrational tran¬ 
sitions, we divide the spectrum into separate regions cov¬ 
ering each OH vibrational transition. This allows indepen¬ 
dent scaling of each transition. We also scale the O 2 emis¬ 
sion around 0.864 /am separately. Furthermore, to account 
for flexure (error in the wavelength calibration) we interpo¬ 
late the spectrum of each transition (OH & O 2 ) forward and 
backwards by one pixel. Ah these spectra are then passed 
to ppxf which finds the best-fit linear summation that re¬ 
produces the skyline residuals while also finding the best fit 
kinematic^] For each galaxy spectrum (whether it be from 
a single bin or a median of the entire field-of-view), we mea¬ 
sure the intrinsic instrument dispersion from the skylines in 
the same wavelength range. 

The SWIFT spectral resolution varies from 40kms -1 
to 65kms -1 . As the resolution of the CvD12 models is 
R=2000 beyond 0.75 /im, it was necessary to correct the 
dispersions found by PPXF for the difference in instrument 
and template resolution: to ah ppxf dispersions we added 
in quadrature the difference between the instrument resolu¬ 
tion in that wavelength range and the template resolution. 


1 We manually altered the internal limits in ppxf to allow nega¬ 
tive sky spectrum weights 
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Figure 2. Plots of the three IMF sensitive spectral features, Nal (left), CaT (centre) and FeH (right) for each pointing in M31 and 
M32. The feature definitions (table [5]) are shown by the vertical dashed lines in each subplot. The spectra have been normalised to the 
continuum bands either side of the feature. M31_l and M32 correspond to median spectra obtained over the entire M31_l and M32 fields. 
For M31_6 in the Nal plot, and for M31_5 and M31_6 in the FeH plot, we show the best fit kinematic template spectra (coloured dashed 
lines) instead of the science data due to contamination from the sky background. Typical ±1 a error bars are denoted by the vertical 
lines at the bottom of each subplot. 


In practice this only significantly affects dispersions less than 
100 km s -1 (M31_6 and M32). 

Using this technique, we are able to clean sky residuals 
from the galaxy spectra (prior to calculating element abun¬ 
dances) using the best-fit ‘second-order’ sky spectrum gen¬ 
erated by ppxf. However, the sky spectrum contains both 
continuum and line emission, and scaling to remove the line 
residuals also scales the continuum. This may be a source of 
bias: adding or subtracting a constant continuum level from 
the galaxy spectrum will alter subsequent equivalent width 
measurements. As there is no reason to believe continuum 
and line emission scale in the same way (unless the con¬ 
tinuum is composed primarily of unresolved faint OH lines 
from the same transition, Davies||2007 ), we subtracted the 
continuum from the second-order sky spectrum before using 
it to remove the skyline residuals in the galaxy spectra. Typ¬ 
ically, the effect of this in equivalent width measurements is 
only a few per cent, but it can be as large as ten per cent. 

For all fields, we fit the kinematics (and second-order 
sky spectra) separately for the Nal-CaT (0.785 /xm to 
0.90 jLtm) and FeH (0.97/xm to 1.02 /xm) regions. For the FeH 
fit, we find it necessary to mask the region between 1.002 /xm 
and 1.008 /xm due to a significant template mismatch be¬ 
tween the CvD12 and M31 spectra. For template spectra, 
we use only the solar metallicity SSPs covering a range of 
ages and IMFs. We do not include spectra for different ele¬ 
ment abundances as these are only available for a 13.5 Gyr 
Chabrier IMF SSP. 

We tried masking the Nal and FeH absorption fea¬ 
tures during the kinematic fitting to test whether fitting 
the second-order sky spectrum affects the strength of these 
absorption features; no significant difference was found for 
data presented here. 


2.5 Index measurements 

Before making index measurements, we de-redshift each 
spectrum to correct for its velocity as determined from the 
kinematic fit, and convolve all our spectra up to a common 
resolution of 200kms _1 using a Gaussian of width given by 


Table 2. Index bandpass and continuum definitions from |Cenarro| 
|et al.| ( |2001| ) and | Conroy &; van Dokkum] ( |2012a| >. Wavelengths 
are in vacuum. We note several different definitions of Nal exist 
in the literature, including Vazdekis et al. (20d2| and|La Barbera| 
|et al.| P 0 T 3 ] ). 


Index 

Feature 

(A) 

Blue Continuum 

(A) 

Red Continuum 

(A) 

NaD 

5878.5-5911.0 

5862.2-5877.2 

5923.7-5949.7 

Nal 

8177.0-8205.0 

8170.0-8177.0 

8205.0-8215.0 

CaT 

8484.0-8513.0 

8474.0-8484.0 

8563.0-8577.0 


8522.0-8562.0 

8474.0-8484.0 

8563.0-8577.0 


8642.0-8682.0 

8619.0-8642.0 

8700.0-8725.0 

PaT 

8461.0-8474.0 

8474.0-8484.0 

8563.0-8577.0 


8577.0-8619.0 

8563.0-8577.0 

8619.0-8642.0 


8730.0-8772.0 

8700.0-8725.0 

8776.0-8792.0 

FeH 

9905.0-9935.0 

9855.0-9880.0 

9940.0-9970.0 


^conv — \/200 2 (cr 2 + cr 2 nstr ), 


( 2 ) 


where cr* in the stellar velocity dispersion measured by ppxf 
and cr instr is the instrumental resolution. Table [2] gives the 
index and continuum definitions for the indices used in this 
paper. 

The spectra around the Nal, CaT and FeH features are 
plotted in Fig. [2] For our resolved data of M31_l and M32, 
we plot a median spectrum from the entire field. We include 
the positions of the spectral index definitions from table [2] in 
each subplot. Fig.[2]shows the high quality of our measured 
spectra. Note that we have omitted the science spectra for 
M31_6 in the Nal plot, and for M31_5 and M31_6 in the FeH 
plot. This is due to strong residual sky lines after optimi¬ 
sation during the sky subtraction stage discussed in section 
2.4. We have instead plotted the best fit kinematic template 
spectra and we show a comparison between the data and 
templates in the appendix. We also note the definition of 
Nal appearing slightly blue-ward of the feature. The feature 
definition is over the Na doublet but there is a TiO band 


next to it, giving the appearance of a larger feature (van 
Dokkum & Co nroy) 2012 ). 
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For M31_1 and M32 we calculate radially binned av¬ 
erage index values and standard errors on the mean from 
the individual spectra in the resolved fields. For M31_2 to 
M31_6 we combine the index uncertainty in quadrature with 
a conservative factor accounting for the uncertainty in tel¬ 
luric correction. 


colour bar indicates the radial distance. Using these index 
maps, we now discuss our measurements in the context of 
our three goals: investigating the disc, bulge and nuclear re¬ 
gion of M31, and the central region of M32. We show model 
variation plots for the VI2 models in the appendix. 


3 RESULTS 


4.1 M31 


Our measurements of the three absorption features as a func¬ 
tion of radius are shown in Fig.[3]and we discuss each feature 
in turn for both galaxies. 

The top panel in Fig. [3] shows our measurements of the 
IMF sensitive index Nal. A clear negative gradient is visible 
in the central ~ 10 x/ of M31, which becomes flat towards the 
edge of the bulge. Our measurement nearest the nucleus is a 
factor of two greater than the average for the bulge. In con¬ 
trast, M32 displays a flat profile of low absorption. It is clear 
from Fig.[3]that, for M31, our best fit template spectra have 
systematically lower values in the central ~ 10" compared 
to the observed spectra. However, in the outer bulge the 
two index measurements are in strong agreement (with the 
exception of M31_6, which we show in the appendix). We 
believe the discrepancy in the central 10 // originates from 
our input template spectra for PPXF not including the SSPs 
with varying element abundances (as these spectra have a 
fixed age and IMF, see section 2.4). 

CaT is shown in the second panel in Fig. [3] M31 displays 
a negative gradient throughout the bulge, dropping from 
7.5 A in the nuclear region down to 6.4 A in the disc. The 
kinematic templates show a small systematic increase of ~ 
0.25 A compared to the measurements within 10 /x , but are in 
very good agreement for the outer bulge measurements. The 
central discrepancy can be partly attributed to differences 
in the continuum measurements between the real and fitted 
spectra - the kinematic fits do not display as strong hydrogen 
absorption in the continuum. Residual telluric absorption 
could also be affecting the continuum level in the observed 
data. M32 displays very high calcium absorption of ~ 8.0 A 
with a slight negative radial gradient. 

Our measurements of FeH are shown in the bottom 
panel of Fig. [3] For both M31 and M32 FeH has a flat pro¬ 
file. The absorption strength is 0.4 A for M31, and between 
0.45 - 0.50 A for M32. 


4 ANALYSIS 

We analyse our measurements of the far red absorption fea¬ 
tures in M31 and M32 using the SSP models of CvD12 and 
V12. Both allow for variations in the IMF slope and age, 
however other parameters are specific to each set. CvD12 
explore a-element enhancement [a/ Fe], as well as variations 
in individual elemental abundance ratios [X/Fe] for 17 dif¬ 
ferent elements (for fixed 13.5 Gyr age and fixed Chabrier 
IMF). V12 explore variations in the total metallicity [Z/H] 
for all IMFs and ages. Note that the VI2 models do not 
cover the FeH feature. Our spectra do not cover the com¬ 
monly used age indicator H/3, so we instead use SDSS g — r 
colour. Fig.[4]shows Nal (left) and FeH (right) against g — r 
for the CvD12 models with the M31 data overplotted. The 


Fig.[4]shows Nal (left) and FeH (right) index measurements 
plotted against SDSS g — r colour as a function of radius. We 
over plot our data onto CvD12 model predictions for both 
indices. Nal appears to suggest an IMF slope around x = 3 
for the central 10" with the bulge values also in a region of 
steepened IMF slope. In contrast, the FeH data follow the 
CvD12 predictions for a Chabrier IMF at solar metallicity 
closely. The discrepancy between Nal and FeH predictions 
can be explained by Na enhancement within the bulge of 
M31. The left plot in Fig.[4]shows that a Chabrier IMF with 
Na enhancement of [Na/Fe] = +0.3 dex could cover the outer 
bulge measurements. The inner bulge and nuclear region 
would require even stronger Na enhancement, with linear 
extrapolation of the models suggesting up to an order of 
magnitude above solar. The g — r data indicate a decreasing 
age between inner and outer bulge along the major axis, 
although the large errors of ±0.03 mag do not allow for any 
stronger conclusions. However, this is in general agreement 
with best fit ages derived by Saglia et al. (2010). 

CaT is known to be inversely proportional to IMF slope 
(see e.g. Cenarro et al. 2003]). However, it is much more 
strongly influenced by [a/Fe] and Ca abundance than by 
IMF. The CvD12 models suggest an order of magnitude in¬ 
crease in the sensitivity of CaT to changes in [a/Fe] com¬ 
pared to changes in IMF (see Fig. A4). The gradient we 
measure is indicative of a steepening IMF slope with ra¬ 
dius. Furthermore, this is in contrast to the FeH measure¬ 
ments. CvD12 note that Na enhancement can change the 
CaT strength. Thus, the flat profile of FeH leads us to con¬ 
clude that the bulge of M31 holds a Chabrier IMF and is 
Na enhanced throughout, more so in the central 10". 

Fig. [3] shows the central 10" of M31 on an index-index 
plot of the two sodium features, Nal and NaD. Model grids 
for CvD12 are fixed at solar metallicity (black), whereas 
those of V12 are at both solar (red) fo r M32 (|Rose et al. 


2005) and [Z/H] = +0.22 (blue) for M31 (Sagl ia et al.|2010| ). 


The radial measurements of the M31 inner bulge and nu¬ 
clear region follow the VI2 model lines of varying age at a 
very bottom-heavy x > 3.3 IMF exponent. Thus the VI2 
models suggest a bottom-heavy IMF in this region, which 
would become Salpeter-like through the rest of the bulge as 
Nal approaches 0.4 A. However, also plotted is a CvD12 SSP 
at 13.5 Gyr, with a Chabrier IMF and sodium enhanced by 
+0.3 dex. It is clear from this SSP that the CvD12 models 
attribute a stronger dependence of NaD on [Na/Fe] than the 
IMF slope, and further it is indicative that the CvD12 mod¬ 
els suggest the level of Na enhancement seen in the central 
10" of M31 is consistent with a normal (between Chabrier 
and Salpeter) IMF. It is not possible to resolve the differ¬ 
ences between the two model predictions without coverage 
of the FeH feature and elemental abundance variations in 
the V12 models. 
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Figure 3. Plots of radial gradients for each index for M31 (black circles) and M32 (blue squares). Shown are the Nal, CaT and FeH 
indices as a function of radial distance. Our measurements of the best fit kinematic template spectra for M31 are shown as the green 
circles. 


4.2 M32 


M32 displays flat profiles for both Nal and FeH, with a 
slightly negative gradient in CaT. Nal is markedly lower 
than in M31 with values around 0.1 A. From Fig [H] both the 
CvD12 and VI2 models suggest this is consistent with a 
Chabrier IMF at solar metallicity around 3-5 Gyr in age. 
The stellar populations of M32 are younger than M31 (see 
O’Connell 1980) and the model predictions for the age 


e.g _ 

from both CvD12 and VI2 are in reasonable agreement with 


long-slit spectroscopy measurements by Rose et al. (2005), 


who measure an age gradient from 3-4 Gyr in the centre to 
6-7 Gyr at 30". A younger stellar population with Chabrier 


IMF is further supported by the FeH index, which is slightly 
stronger than in M31. 

A striking difference is the much higher CaT values 
compared to M31. This has previously been noted by Ko- 


rmendy & Bender (1999) and Cenarro, Cardiel V Gorgas 


(2008). Thomas, Maraston V Bender (2003) suggest this 


can be explained by the fact that M32 is younger than 
M31, and the younger stars contaminate the CaT index, 
which is known to be affected by Paschen lines. However, 
the CaT* index is defined by Cenarro et al. (2001) as CaT* 


= CaT — 0.93 x PaT, where PaT is an index defined over 
several Paschen lines next to CaT, in order to calibrate out 
the hydrogen lines. We measured the CaT* index within the 
central 2" of M32 as CaT* = 7.39 zb 0.22 A and within the 
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NaI0.82 [A] FeH0.99 [A] 


Figure 4. Plots of Nal (left) and FeH (right) against SDSS g — r colour showing CvD12 models with M31 data overplotted. B-L 
corresponds to a bottom-light IMF from van Dokkum (2 008| ), and Chab to a[Chabrier (2 Q03| > IMF. The colour bar corresponds to the 
radial distance of each pointing to the centre of M31. Symbol sizes also decrease with increasing radial distance. Colour error bars are 
omitted from the points for clarity and are shown in the bottom right of each plot. 


same aperture of M31 as CaT* = 6.75 d= 0.16 A, which still 
shows greater absorption for M32. It would be an interest¬ 
ing test to see whether other dwarf galaxies of the same age 
as M32 also display strong calcium absorption. Overall our 
measurements of M32 are fully consistent with a Chabrier 
IMF. 


5 DISCUSSION 


The form of the IMF in the nuclear region of M31 has been 
of intense historical debate. The main factor driving the idea 
of a bottom heavy IMF has been the prominent Nal feature 
in the central ~ 3". Spinrad & Taylor (1971) interpreted 


this as the presence of a large number of dwarf stars. Fur¬ 
ther observations by many groups (see e.g. Whitford 1977 


Fro gel et al.|1978[ Cohen|1 978| Faber & French 1980; CajTer^ 

Visvanathan &; Pickles||1986| |Alloin &; Bica||1989| |Delisle 


Hardy 19 92|) hav e led to conflicting results. [Conroy Sz van 

Dokkum (2012b) observed the central 15 pc of M31 (4 ) 
and inferred a mass-to-light ratio consistent with an IMF 
slope between Chabrier and Salpeter. They concluded that 
the central region is entirely consistent with a Milky Way 
IMF, albeit requiring very large sodium enhancement of the 
order [Na/Fe]~ 1.0. 

All current work suggesting increasing IMF slopes with 
velocity dispersion in the centres of galaxies has been for 
massive early-types, which gene rally have larg e dispersions, 
upwards of 300kms _ 


Ferrer as et al. 


(2013) 


La Barbera 


et al. (2013) and Spiniello et al. (2014) find similar em¬ 
pirical relations between the IMF slope and velocity dis¬ 
persion. M31 has a dispersion of a ~ 160kms _1 (out to 
0.4kpc which we adopt as a typical central dispersion), 
which gives an IMF slope of x ~ 1.9 from all three rela¬ 
tions and is consistent with a normal (i.e. not bottom-heavy) 


IMF. We measure a Chabrier IMF throughout M31 even at 
the largest radius which has a dispersion of cr ~ 90kms _1 , 
and also in the centre of M32 which has a central dispersion 
of ~ 60kms _1 ; this suggests that the relations linking IMF 
to velocity dispersi on shou ld become f i xed at Chabrie r IMF 
for cr < 150kms _1 . Chabrier, Hennebelle & Chariot (2014) 


have presented a theoretical scenario for dwarf enrichment 
and a steepening of the IMF within the framework of grav¬ 
itationally collapsing turbulent molecular clouds. They de¬ 
duce an upper limit to the IMF slope of x ~ 2.8 in extremely 
high density and turbulent systems. We note that the cen¬ 
tres of M31 and M32 have very high stellar density and yet 
we measure a Chabrier IMF. 

The physical size scales here are also worth consider¬ 
ing. The central 8" of M31 corresponds to 30 pc, which is a 
very small region of space in which to create and maintain 
an overabundance of dwarf stars (or in fact, a population 
with [Na/Fe]~ +1.0). Furthermore, sodium excess in stars 
(Bril ey et al.|1 996) and in galaxies ( |Jeong et al.|2013 ) is not 
particularly rare, suggesting that M31 may not be peculiar. 
In fact, Jeong et al. (2013) measure NaD absorption as high 


as ~ 7 A in some ETGs and late type galaxies (LTGs); much 
higher than seen in M31. They conclude that the NaD in¬ 
dex in most Na enhanced LTGs is strongly contaminated by 
the interstellar medium (ISM) and/or dust. These findings 
caution against using sodium as a sole IMF indicator, es¬ 
pecially for late type galaxies. Most importantly, our mea¬ 
surements of the FeH absorption reject the possibility of 
a dwarf -dominated nuclear region. |Conroy &; van Dokkum| 
(2012b) also reached this conclusion for the central 4" and 


our results have extended that conclusion throughout the 
bulge and into the disc. 

It is an important point that the two SSP models give 
different predictions when considering the sodium indices as 


highlighted by Spiniello, Trager Sz Koopmans (2015). Specif- 
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Figure 5. SSP model grids for CvD12 (black) and V12 (red and blue) showing variations in Nal against NaD. The CvD12 models are at 
solar metallicity, whereas the V12 models are at both solar (red, for M32) and metal enhanced with [Z/H]=+0.22 (blue, for M31). Note 
that the varying age vectors from the VI2 models are parallel to the [Na/Fe] vector from the CvD12 models (black arrow). The central 
10 // of M31 (circles) and M32 (squares) are plotted. Symbol sizes decrease with increasing radial distance. We use the NaD measurements 
of |Davidge| (1991 1 997} that cover r ^ 10". IMF labels are the same as in Fig. [ 4 ] 


ically VI2 attribute a much stronger dependence of NaD 
on the IMF slope than the CvD12 models. This is clear in 
Fig. [5] and gives apparently discrepant predictions for the 
IMF slope in M31. We conclude with |Spiniello, Trager <V| 


Koopmans (2015) that sodium is complex and poorly un¬ 


derstood. However, the V12 models do not cover the FeH 
wavelength range so no constraint on the IMF can be made 
using this feature with the VI2 models. 


6 CONCLUSIONS 

Using the Oxford SWIFT integral field spectrograph we have 
undertaken a study of M31 and M32 with the aim of discern¬ 
ing possible gradients in IMF-sensitive far red absorption 
features. We obtained high S/N integral field spectroscopy 
of M31, covering the nuclear region, bulge and old disc at 
~ 700 // (2.7kpc), as well as for the central 10"(38 pc) of M32. 
We conclude the following: 

1. In M31, the Nal index shows a flat equivalent width 
profile at ~ 0.4 A through the majority of the bulge, with a 
strong gradient up to 0.8 A in the central 10". FeH is mea¬ 
sured to be flat at 0.4 A on all scales. CaT shows a gradual 
positive gradient inwards through the bulge. M32 displays 
flat profiles for all three indices: slightly stronger FeH ab¬ 
sorption than M31, anomalously strong CaT absorption in 
agreement with other studies, and very weak Nal absorp¬ 
tion. 

2. M31 is consistent with a Chabrier IMF throughout 
the central region, bulge and old disc from analysis of the 
hat FeH index using the CvD12 models. We also infer a 


gradient in sodium enhancement of [Na/Fe] ~ +0.3 dex in 
the disc and outer bulge, rising up to [Na/Fe] ~ +1.0 dex in 
the nuclear region. 

3. Our measurements of Nal and FeH for M32 are con¬ 
sistent with a Chabrier IMF for ages ~3-5Gyr using the 
CvD12 and VI2 models. 

4. The use of sodium indices alone is a poor method for 
investigating IMF variations as it is susceptible to sodium 
enhancement and [a/Fe] variations. Sodium is a poorly un¬ 
derstood element as highlighted by the discrepant predic¬ 
tions of the CvD12 and V12 models. 

5. The centres of M31 and M32 have very high stel¬ 
lar densities and yet we measure a Chabrier IMF in both 
regions. 
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APPENDIX A: ADDITIONAL PLOTS 

We show here several additional plots concerning the ac¬ 
curacy and limitations of our telluric correction and sky 
subtraction techniques, as well as comparison plots between 
CvD12 and VI2 models. 

Fig. |Al| shows examples of our telluric correction around 
the Nal feature for M31_2 and M31_6. In each plot we use our 
kinematic template fits as an assumed ‘telluric free’ galaxy 
spectrum. The residuals after telluric correction are 0.8 per 
cent and 1 per cent rms respectively. We also check this 
for CaT and FeH and find residuals < 1 per cent for these 
features. 

Fig. |A2| shows the M31_6 spectrum around the Nal fea¬ 
ture. The second-order sky subtraction procedure has deep¬ 
ened the feature with extra sodium absorption, leading to 
a deeper feature. We checked this for all spectra to ensure 
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Figure Al. Telluric correction for M31_2 (left) and M31_6 (right). The black line shows the ratio between original spectrum and 
kinematic template (assuming this is the intrinsic ‘telluric free’ spectrum), the red line shows the same ratio with telluric division, and 
the green line shows the telluric spectrum scaled by factor of 0.5 and plotted at 1.3 for clarity. 



A [jum] 


Figure A2. Plot showing the M31_6 spectrum before (black) 
and after (red) the 2nd order sky subtraction routine (section 
2.4). The procedure has deepened the Nal feature with (false) 
excess sodium emission, which leads to a larger index value from 
the convolved spectra (dashed lines). 


systematic errors were minimised and only found this con¬ 
tamination in the M31_6 spectrum. We show this plot to 
highlight why we have used the best fit kinematic template 
to measure Nal for M31_6 instead of the science spectrum. 

Fig. |A3| shows our M31_l, M31_5 and M31_6 science 
spectra around the FeH features with their best fit kinematic 
templates. The presence of residual sky lines after the second 
order sky subtraction stage in the M31_5 and M31_6 spectra 
leads us to use the template spectra to measure the FeH 
feature for these fields instead of the science spectra. 

Fig. |A4| shows comparison model grids between CvD12 
(left) and V12 (right) for Nal (top row), CaT (second row) 
and FeH (bottom; CvD12 only) features. We plot each index 
against variation in SDSS g — r colour. These maps serve to 
show and compare the behaviour of each set of SSP mod¬ 
els for each index within their respective model parameter 
spaces. Both sets of models qualitatively agree in their be¬ 
haviour for variations in IMF and age. However quantitative 
differences exist, especially for the CaT index. 




Figure A3. Plotted spectra showing science spectra (black lines) 
with best fit models (red lines) for three fields, M31_l, M31_5 
and M31_6. Solid lines show the spectra at their intrinsic res¬ 
olution, and dashed lines show the spectra convolved up to 
200kms —1 . The spectra have each been normalised over the 
pseudo-continuum of the FeH feature. 


This paper has been typeset from a TgX/ BT^X file prepared 
by the author. 
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Figure A4. Plots of Nal index (top row) and CaT (second row) against SDSS g — r colour showing model variation for CvD12 (left) 
and V12 (right) models. Also shown is FeH (bottom) against g — r for CvD12 models only, as FeH is not covered by V12. For the 
V12 models, black lines and circles are different IMFs at 14 Gyrs and [Z/H]=0.0, red squares are different ages, and blue diamonds are 
different metallicities. For the CvD12 models, black lines and circles show different IMFs, red lines and squares show different ages, blue 
lines and diamonds show varying [a/Fe] abundance, and green dashed lines and diamonds show varying [X/Fe]. 


















